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Tau-pair process has been analyzed in the ILD detector model as a benchmark process 
for Lol. Results of background rejection, forward-backward asymmetry and polarization 
measurements are obtained with full detector simulation. Impact of detector parame- 
ters for tau-pair analysis is also discussed in this paper. 



1 Motivations of tau-pair study 

Tau-pair process (e+e^ ^Z* t+t^) at y/s = 500 GeV is one of the benchmark processes^ 
proposed by Research Director. According to the report, this process is a good sample to 
examine detector performances of 

• tau reconstruction, aspects of particle flow, 

• ttq reconstruction, 

• tracking of very close- by tracks. 

In this process, tau leptons are highly boosted (7 ^ 140), thus decay daughters (mainly 
charged and neutral pions, muons and electrons) are concentrated in a very narrow angle. 
Reconstruction of ttq from two photons is especially challenging for the ILC detectors, and 
much depends on detector parameters, so it is a good measure for detector optimization. 

Required observables arc cross section, forward-backward asymmetry and polarization of 
tau leptons. The polarization measurement requires identification of tau decays, including 
reconstruction of ttq. Efficiency and purity of event selection cuts should also be used for 
comparison of detector performances. 

For physics motivation, tau-pair process is important as a precision measurement of the 
electroweak theory. For example, measuring cross section and forward-backward asymmetry 
of tau-pair process very precisely can probe existence of heavy Z' boson. 
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Geometry 


gldapr08_14m 


gldprim_v04 


j41dc_v04 


LDCPrime.02Sc 


Software 


Jupiter 


Jupiter 


Jupiter 


Mokka 


Magnetic field 


3 Tesla 


3.5 Tesla 


4 Tesla 


3.5 Tesla 


TPC Rmin 


43.7 cm 


43.5 cm 


34.0 cm 


37.1 cm 


ECAL Rmin 


210 cm 


185 cm 


160 cm 


182.5 cm 


ECAL thickness 


19.8 cm 


19.8 cm 


19.8 cm 


17.2 cm 


HCAL thickness 


120 cm 


109 cm 


96 cm 


127.2 cm 


ECAL granularity 


1x1 cm 


1x1 cm 


1x1 cm 


0.5x0.5 cm 



Table 1: Detector geometries used in this study. 



2 Analysis framework and events 

2.1 Monte Carlo simulation and detector geometries 

The ILD group has two full detector simulation models, Mokka and Jupiter. Mokka origi- 
nates in LDC detector and Jupiter originates in GLD detector, and both are based on Geant4 
Monte Carlo simulation. In this study I used both simulation models. Mokka has geometries 
with detailed implementation of detector components. I used simulated events processed in 
Mokka LDCPrime_02Sc geometry. The ILD group has simulated quite a large fraction of 
full Standard Model (SM) samples, required in the benchmark report, in LDCPrime_02Sc 
geometry and I used the events to estimate and optimize background suppression. In con- 
trast, Jupiter has relatively rough geometries and full SM samples have not been processed, 
but we have tau-pair samples in several Jupiter geometries with three detector parame- 
ters, gldapr08_14m, gldprim_v04, j41dc_v04. Summary of detector geometries is shown in 
Fig. [ij In rough summary, three geometries in Jupiter differ in sizes, gldapr08_14m is large, 
gldprim_v04 is middle, j41dc_v04 is small in size. Magnetic field is such that BR for each 
geometry is almost the same. LDCPrimc_02Sc is almost as same as gldprim_v04 in size, but 
it has a finer ECAL granularity of 0.5x0.5 cm. Detailed geometry is much different between 
Jupiter and Mokka geometries. 

For event reconstruction, including smearing of tracker and calorimeter hits, tracking, 
clustering and particle flow, I used MarlinReco framework with PandoraPFA particle flow 
algorithm. Raw output of Jupiter is not compatible with MarlinReco LCIO format, but we 
have a converter to obtain LCIO files of Jupiter events. PandoraPFA is especially optimized 
for Mokka geometries, so particle flow performance of Jupiter is slightly worse. 

2.2 Event samples 

For signal tau-pair events, we use events generated in DESY. Whizard 1.51 and TAUOLA 
are used to generate the events. SLAG standard samples for Lol are not used because of 
polarization issues. I used 80 fb~^ signal sample of each geometry for ApB and polarization 
analysis without background. For background study, SLAG standard samples are used. All 
events simulated in LDCPrimc_02Sc geometry, about 20 million events in total, are processed 
with my analysis cuts. 

Integrated luminosity is assumed to be 500 fb"-'^ each for two polarization setups, ej^ej^ 
and epC^. Assumed polarization ratio is 80% for electron and 30% for positron (i.e. for 
^L^R setup 90% of electrons arc Icftly polarized and 65% of positrons arc rightly polarized). 
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2.3 Tau clustering 



For tau clustering, an original clustering processor (TaJet) is applied to the output of Pan- 
doraPFA. Following is a procedure of the processor. 

1. Sort particles in energy order. 

2. Select the most energetic charged particle (a tau candidate). 

3. Search particles to be associated to the tau candidate. Criteria is: 

(a) Opening angle to the tau candidate is smaller than 50 mrad.. or 

(b) Opening angle to the tau candidate is not larger than 1 rad. and invariant mass 
with the tau candidate is less than 2 GcV (nv = 1.777 GcV). 

4. Combine energy and momentum of the tau candidate and associated particle and treat 
the combined particle as the new tau candidate. 

5. Repeat from 3. 

6. After all remaining particles do not meet the criteria, remaining most energetic charged 
particle is the next tau candidate. (Repeat from 2.) 

7. After all charged particles are associated to tau candidates, remaining neutral particles 
arc independently included in the cluster list as neutral fragments. 

In the clustering stage, events with > 6 tracks arc pre-cut to accelerate clustering since 
> 99% of tau decays have < 3 charged particles. Event with only one positive and one 
negative tau clusters are processed with latter analysis. 



3 Background suppression 

Main background of tau-pair analysis is Bhabha (c+c^ — s-e+e"), WW — > ivlv and 77 
T~^T~ . Since cross sections of Bhabha and two-photon events are huge (about 10^ and 10^ 
larger than signal, respectively), we need tight selection cuts for those background events. 
Following cuts are applied to signals and all SM background events after the tau clustering. 

1. Number of tracks < 6. Included as a pre-cut in tau clustering processor. 

2. Only one positive and one negative tau clusters must exist in the event. 

3. Opening angle of two tau candidates must be > 178 deg. 
This cut efficiently suppresses WW ivlv background. 

4. ee and /i^ events are rejected. 

Charged particles depositing > 90% of their energy in ECAL are identified as electrons, 
and charged particles depositing < 70% of their energy (estimated by curvature of their 
tracks) in ECAL+HCAL are identified as muons. Events with two electrons or two 
muons are rejected in this cut. This cut is especially for suppressing Bhabha and 
e"'"e~ — > Ijl^ I-l" events. Signal loss is about 6%. 
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5. I cos6'| < 0.9 for both tau clusters. 

t-channel Bhabha events are almost completely suppressed by this cut. 20% of signal 
events are lost. 

6. 40 < Evis < 450 GeV. 

Lower bound suppresses 77 t^t^ events, and upper bound suppresses Bhabha 
events. Signal lost is negligibly small. 



Table [2] shows the result of 
these cuts. 77 t^t~ back- 
ground is currently not included, 
but generator-level study shows 
the effect of the 77 — > r+r" back- 
ground is not significant (most 
events are eliminated by the Evis 
cut). Number of events are nor- 
malized to 500 fb~^. ej^ej^ po- 
larization (80% and 30%, respec- 
tively) is assumed. In total, num- 
ber of background events is about 10% of signal, not significant. If we assume that we know 
the shape of the background, statistical error of the background is much smaller than signal 
statistics and negligible for the further study. 77 t^t^ background is planned to be 
included in the Lol study. 



Cuts 


Signal 


Background 


# tracks, # clusters 


5.7 X 10^ 


7.9 X 10** 


Opening angle > 178 deg. 


1.6 X 10^ 


1.3 X 10** 


1 cos 6*1 < 0.9 


1.3 X 10^ 


1.2 X 10^ 


ee, /ijU veto 


1.2 X 10^ 


6.2 X 10^ 


45 < Evis < 450 


1.2 X 10^ 


1.3 X 10"^ 



Table 2: Cut statistics for background suppression. 



4 Forward-backward asymmetry 



Figure [T] shows a simulated result on an- 
gular distribution of r"*" leptons in LD- 
CPrime_02Sc model. All events passed cri- 
teria described in the previous section are 
put into the histograms. Clear asymmetry 
can be seen in signal distribution. There are 
several bins where background is very large, 
but this is a result of low MC statistics of 
Bhabha, 7*7* and 67* processes (about 0.1 
fb~^ each). We plan to improve MC statis- 
tics of those events by applying preselec- 
tions at generator level. 

Forward-backward asymmetry can 
calculated by 
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Figure 1: Angular distribution of r"*" momen- 
tum direction. Black crosses show signal and 
red crosses show background. Vertical axis is 
normalized to 500 fb~^. Error bars represent 
errors in current MC statistics. 
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where Nb is number of events in back- 
ward region {cos 6 < 0) and Np is num- 
ber of events in forward region (cos 6 > 0). By estimated number of signal events 
(Np = 8956, TVs = 2893), Afb is estimated to be 51.2±0.25% (background statistics is 
not included). Since total estimated number of background events is about 10% of sig- 
nal, effect of background to statistical error is smaller than signal statistics if background 
distribution can be well determined. 

5 Polarization analysis 
5.1 Event selection 

There are five dominant decay modes of tau leptons, t"*" — > e'^V^VT (17.9%), r"*" — > ^l^IT^Vt 

(17.4%), r+ TT+Vr (10.9%), T+ -> p+Vr ^ TT+Tr°Vr (25.2%), and T+ -> atvr TT+Tr°Tr°Vr 

(9.3%). 

Since first two modes are leptonic 3-body decay and polarization information is partially 
lost due to the missing neutrinos, and the last ai mode has relatively low branching ratio, 
we currently use only r+ -k^Vt and r+ p^Vr modes. These modes are selected by 
following criteria. 

1. Tau clusters with one charged tracks are selected. 

2. Clusters with electrons and muons are eliminated. Muon identification is the same 
as that in SM suppression cut. Electron identification is ECAL deposit > 90% for 
r+ n^Vr selection and > 97% for r+ p'^Vr selection. 

3. Clustered with energy > 10 GcV is eliminated (since lepton identification is currently 
poor in low energy clusters). 

4. Check whether neutral particles are associated in the cluster. If > 1 GeV neutral 
particles are not associated, the cluster is treated as a t+ ^ n^Vr event. If > 10 GeV 
neutral particles are associated, the cluster is treated as a r+ ^ p^i'r event candidate. 

5. For P^Vt candidates, invariant mass of p is calculated from 4-momenta of 
charged pion and whole cluster. Clusters with invariant mass around 200 McV to nip 
(570 to 970 MeV) are accepted. 

6. Optional tt^ mass cut is applied to the r+ p^Vr candidates. In this cut, invariant 
mass of tt" is calculated with clusters which have > 2 neutral particles Events with 
invariant mass of to 200 MeV (m° = 135 MeV) are accepted. All events with only 
one neutral particle are eliminated by this cut (if applied). 

Figure [2] and [3] shows invariant mass distribution of p and tt*' with gldapr08_14m (noted 
GLD in the plot), gldprim_v04 (GLD'), j41dc.v04 (J4LDC) and LDCPrime.02Sc (LDC) 
geometries. Especially tt*^ invariant mass distribution shows large difference between ge- 
ometries. Larger and higher granularity geometry apparently gives better results in 7r° 
reconstruction. 

Table [3] shows a result of selection efficiency and purity for each detector geometry. For 
r"*" — > p'^Vr mode, both results of selection without tt'' invariant mass cut and selection with 
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Figure 2: Invariant mass distribution of p 
with four geometries. 
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Figure 3: Invariant mass distribution of n' 
with four geometries. 



Geometry 


gldapr08_14m 


gldprim_v04 


j41dc_v04 


LDCPrime.02Sc 


TT mode efficiency 


21.3% 


21.4% 


21.4% 


21.2% 


TT mode purity 


85.7% 


83.6% 


80.8% 


88.5% 


p mode eff. wo/ tt" cut 


12.7% 


12.1% 


11.3% 


12.8% 


p mode pur. wo/ tt" cut 


83.4% 


81.8% 


81.4% 


85.7% 


p mode eff. w/ tt*^ cut 


5.31% 


4.32% 


3.72% 


6.38% 


p mode pur. w/ tt'^ cut 


92.3% 


90.3% 


90.5% 


93.9% 



Table 3: Selection efficiency and purity for polarization analysis. 
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Figure 4: Angular distribution of charged pion in tau-momentum frame, left: e^ejlj^ polar- 
ization, right: e^e^ polarization. 
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Figure 5: Distribution of parameter y. See text for details of y. left: mode selection without 
TT^ invariant mass cut, right: selection with tt'^ mass cut. 



7r° invariant mass cut arc shown in parallel, tt'^ mass cut gives higher purity in event selection 
but efBciency becomes much less, thus the cut seems not practical for analysis in current 
geometries. For all selection LDCPrimc_02Sc gives the best result, and gldapr08_14m follows 
the next. 

5.2 Polarization measurement 

To obtain polarization of t+ tt^Vt events, cos0 distribution of pion momentum direction 
with respect to tau momentum direction should be observed. 

Figure[4]shows the cos^ distribution. For the e^ Cj!;^, number of events is larger in cos 9 < Q 
area, and for the epC^, number of events is larger in cosO > area. Polarization of tau 
Icptons can be determined by the ratio of number of events between left and right half of the 
graph, or by linear fit of the histograms. Analysis shows the polarization can be determined 
by 1.2-1.3% statistical error, but the amount of remaining background varies by geometries 
as a result of difference in selection purity (for detailed numbers, see the slide[T]). 

For analysis of t+ — > p^v-r polarization, analysis is more complicated. To obtain po- 
larization of tau leptons, we can use p polarization, indicated by angular distribution of 
p^ <— Tr+TT*^ decay, in addition to p angular distribution with respect to t frame. To com- 
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bine those indicator, we use y parameter, defined in [3] as 



According to [3] , polarization of tau leptons can be determined by 

Figure [5] shows the y distribution. From the distribution, can be determined by 1.1- 
1.2% statistical error if we do not apply tt^ mass cut, and 1.7-2.3% with 7r° mass cut (again, 
detailed number can be seen in the slide[l]). Difference can be observed between geometries, 
due to selection efficiency, but with no mass cut the difference is not so large. 



6 Summary 

Tau-pair process has been analyzed in the ILD detector models. It is found that tau-pair 
forward-backward asymmetry observation of better than 1% resolution can be achieved in 
the current ILD detector models and no large difference between detector models are found. 
For the polarization analysis, clear dependence is seen in ttq reconstruction. LDCPrime_02Sc 
gives the best result and gldapr08_14m follows. With the result we can estimate that larger 
and more granular detector models give better results. 
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